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Hie indentation load-displacement behavior of six materials tested with a Berkovicb 
indenter has been carefully documented to establish an improved method for determining 
hardness and clastic modulus from indentation load-displacement data. The materials 
included feed silica, soda-lime glass, and single crystals of aluminum, tungsten, quartz, 
and sapphire. It is shown that the load-displacement curves daring unloading in these 
materials are not linear, even in the initial stages, thereby suggesting thai the flat punch 
approximation used so often in the analysis of unloading data is not entirely adequate. 
An analysis technique is presented that accounts for the curvature in the unloading data 
and provides a physically justifiable procedure for detennining the depth which should 
be used in conjunction with the indenter shape function to establish the contact area at 
peak load. Hie hardnesses and clastic moduli of the six materials are computed using 
the analysis procedure and compared with values determined by independent means to 
assess the accuracy of the method. The results show that with good technique, moduli 
can be measured to within 5%. 



I. INTRODUCTION 

Great strides have been made over the past few 
years in the development of techniques for probing the 
mechanical properties of materials on the subrnicron 
scale. 1 ' 3 The advances have been made possible by the 
development of instruments that continuously measure 
force and displacement as an indentation is made. 4 " 7 The 
indentation load-displacement data thus derived can be 
used to determine mechanical properties even when the 
indentations are too small to be imaged conveniently. 
Since the indentation positioning Capability of some of 
the instruments h also in the subrnicron regime, a means 
is available by which the mechanical properties of a 
surface can be mapped with subrnicron resolution. An 
instrument with such capabilities can be described as a 
mechanical properties microprobe^ 10 

The two mechanical properties measured most fre- 
quently using load and depth sensing indentation tech- 
niques are the elastic modulus, E 7 and die hardness, 
H. In a commonly used method, data are obtained 
from one complete cycle of loading and unloading.* 
The unloading data are then analyzed according, to a 
model for the deformation of an elastic half space by 
an elastic punch which relates the contact area at peak 
load to the elastic modulus. Methods for independently 
estimating the contact area from the indenter shape func- 
tion are then used to provide separate measurements of 
E and A. 



The elastic contact problem, which plays a key 
role in the analysis procedure, was originally considered 
in the rate 19th century by Boussinesq u and Hertz. 12 
Boussinesq developed a method based on potential the- 
ory for computing the stresses and displacements in an 
clastic body loaded by a rigid, axisymmetric indenter. 
His method has subsequently been used to derive so- 
lutions for a number of important geometries such as 
cylindrical and conical indenters. 13 *' 4 Hertz analyzed the 
problem of the elastic contact between two spherical sur- 
faces with different radii and elastic constants. His now- 
classic solutions form the basis of much experimental 
and theoretical work in the field of contact mechanics 15 
and provide a framework by which the effects of non- 
rigid indenters can be included in the analysis. Another 
major contribution was made by Sneddon, who derived 
general relationships among the load, displacement, and 
contact area for any punch that can be described as a 
solid of revolution of a smooth function. 14 ' 17 His results, 
which will be used extensively in this paper, show 
that the load-displacement relationships for many simple 
punch geometries can conveniently be written $$ 



P = *h m 



(1) 



where P is the indenter load, h is the elastic displacement 
of the indenter, and a and m are constants. Values of 
the exponent m for some common punch geometries are 
m = 1 for flat cylinders, m = 2 for cones, m =s 1.5 tor 
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spheres in the limit of small displacements, and m = 1.5 
for paraboloids of revolution. 

Modeling indentation contact in a way (hat includes 
plasticity is a much mote complex problem. Since the 
constitutive equations are nonlinear and a number of ma- 
terial parameters must be included to describe material 
behavior (e.g^"yield strength and work hardening coeffi- 
cient), analytical solutions are not easily obtained. 15 As 
a result, much of our undemanding of the importance of 
plasticity in indenter contact problems has been derived 
through experimentation and finite element simulation. 

The earliest experiments in which load and displace- 
ment sensing indentation methods weir used to mea- 
sure mechanical properties were performed by Tabor/ 8 
who studied the indentation of a number of metals 
deformed by hardened spherical indenters. A similar 
Study was subsequently undertaken by Still well and 
Ibbor to examine Che behavior of conical indenters,' 9 
One particularly important observation resulting from 
these studies concerns the shape of the hardness im- 
pression after the indenter is unloaded and the material 
elastically recovers. The experiments revealed that, at 
least in metals, the impression formed by a spherical 
indenter is still spherical with a slightly larger radius 
than the indenter, and the impression fonned by a 
conical indenter is Still conical with a larger included 
tip angle. The importance of these experiments is that 
since elastic contact solutions exist tor each of these 
geometries (i.e., a spherical indenter in a spherical hole 
and a conical indenter in a conical hole), the ways 
in which plasticity affects the interpretation of elastic 
unloading data can be dealt with by taking into account 
the shape of the perturbed surface in the elastic analysis. 
Tabor used these results to show that the shape of the 
entire unloading curve and the total amount of recovered 
displacement can be accurately related to the elastic 
modulus and the size of the contact impression for 
both spherical and conical indenters. Other important 
observations resulting from these studies include (1) 
the diameter of the contact impression in the surface 
formed by conical indenters does not recover during 
unloading— only die depth recovers; (2) the indentation 
must be loaded and unloaded a few times before the load- 
displacement behavior becomes perfectly reversible; i.e., 
a limited amount of plasticity sometimes occurs in each 
of the first few loadings and unloadings; and (3) effects 
of non-rigid indenters on the load-displacement behavior 
can be effectively accounted for by defining a reduced 
modulus. Er, through the equation 

where E and u are Young's modulus and Poisson's ratio 
for the specimen and E t and Vi are the same parameters 
for the indenter. 



Interest in load and displacement sensing indentation 
testing as an experimental tool for measuring elastic 
modulus began in the early 1970*s with the work of 
Bulychev, Alekhin, Shorshorov, and co-workers. 20 - 24 
These investigators used instrumented mkrohardncs* 
testing machines to obtain indentation load-displacement 
data tike that shown schematically in Fig. i which was 
then analyzed according to the equation 

Here, S = dP/dh is the experimentally measured 
stiffness of the upper portion of die unloading data, E r 
is the reduced modulus (previously defined), and A is 
the projected area of the clastic contact. By measuring 
the initial unloading stiffness and assuming that the 
contact area is equal to the optically measured area of the 
hardness impression, the modulus can thus be derived. 

Equation (3) has its origins in elastic contact the- 
ory. While originally derived for a conical Indenter, 
Buytchcv & aL showed that Eq. (3) holds equally well 
for spherical and cylindrical indenters and speculated 
that Eq. (3) may apply to other geometries as well. 21 
Pharr, Oliver, and Brotzcn have subsequently shown that 
Eq. (3) applies to any indenter that can be described as 
a body of revolution of a smooth (unction. 25 Bulychev 
ei at, also argued that significant deviations from the 
behavior predicted by Eq, (3) Should not occur for 




DISPLACEMENT, h 

FIG. 1. A schematic representation of load versus indemcr displace- 
ment data for an indentation experiment The quantities shown are 
P mtM : the peak indentation load; A***: the indenter displacement at 
peak load; «y; the ftoat depth of the contact impression after unloading; 
and S: the initial Uploading *liffnctt. 
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pyramidal indented. TTxey noted good agreement be- 
tween generally accepted values of moduli and values 
computed from Eq. (3) using indentation data obtained 
with a Vickers indenter. That Eq. (3) works well for at 
least some mdenters that cannot be described as bodies 
Of revolution has been confirmed by the finite element 
calculations of King.** These calculations show that the 
deviations from Eq. (3) for fiat ended punches with 
square and triangular cross sections are only 1-2% and 
3.4%, respectively. 

In the early 1980'$, it was realized that load and 
depth sensing indentation methods could be very useful 
in the measurement of die mechanical properties of 
very thin films and surface layers, and instruments for 
producing suhmicron indentations were developed.*" 7 
For practical reasons, some means other than direct 
observation of the hardness impressions was needed 
to measure contact areas since imaging very small in- 
dentations is both time-consuming and difficult Oliver, 
Hutchings, and Pethica suggested a simple method based 
on measured indentation load-displacement curves and 
a knowledge of the indenter area function (or shape 
function), that is, die cross-sectional area of the indenter 
as a function of the distance from its tip. 1 * 2 The method 
is based on the notion that, at peak load, the material 
conforms to the shape of die indenter to some depth; if 
this depth can be established from the load-displacement 
data, the projected area of contact can be estimated 
directly from the shape function, 'two obvious choices 
for the depth, though not the only ones, are the depth 
at peak load, /w (Le., the maximum displacement in 
the loading cycle), and the final depth, h f (Le. t die 
residual depth of the hardness impression after final 
unloading), both of which are easily determined from 
indentation load-displacement data (see Fig. 1). Using 
TEM replication methods to establish the shape func- 
tion, Oliver er al. found that the final depth gives a 
better estimate of the contact area than the depth at 
peak load. 2 

Doemer and Nix subsequently put together many of 
these ideas to produce the most comprehensive method 
to date for determining hardness and modulus from 
indentation load-displacement data. 3 Their approach is 
based on the observation that during the initial stages 
of unloading, the elastic behavior of the indentation 
contact is similar to that of a flat cylindrical punch; 
that is, the area of contact remains constant as the 
indenter is unloaded. This was justified by experimental 
observations which suggest that, for some materials, 
the initial portions of unloading curves arc linear, as 
would be expected for the flat-punch geometry. The 
unloading stiffness dP/dh is then related to the modulus 
and contact area through a relationship equivalent to 
Eq. (3). To evaluate independently the contact area, 
Doemer and Nix proposed a simple empirical method 



based on extrapolating the initial linear portion of the 
unloading curve to 2ero load and using the extrapolated 
depth with the indenter shape function to determine 
the contact area. Careful experiments in one material, 
METGLAS® 2826, confirmed that the extrapolated depth 
gives a better estimate of the contact area than either the 
depth at peak load or the final depth. This observation 
was later confirmed with finite element simulations of the 
indentation of silicon and nickel by corneal indenters. 27 
With the contact area so determined, die modulus can be 
computed from Eq. (3) and the hardness from its normal 
definition: 



H = 



where P m h the peak indentation load and A is the 
projected area of the hardness impression. As a practical 
matter, Doemer and Nix suggested that die unloading 
stiffness can be computed from a linear fit of die upper 
one-third of the unloading curve. 

One disturbing feature of the Doerner-NIx method 
concerns their assumption of linear unloading. Dur- 
ing the past several years, we have conducted load 
and displacement sensing indentation tests on a large 
number of materials, and out observations have led 
us to believe that unloading curves are rarely, if ever, 
linear, even in the initial stages of unloading. Rather, 
unloading data are better described by power laws like 
Eq. (1) with exponents ranging from about 12 to 1.6. 
In addition, by employing a special dynamic technique 
by which stiffness can be measured continuously dur- 
ing indentation, we have found that unloading con- 
tact stiffnesses change immediately and continuously 
as the indenter is withdrawn, as would be expected 
from continuous changes in contact area. Together, these 
observations suggest that the flat punch approximation 
is not an entirely adequate description of real material 
behavior. 

In this paper, a new method for analyzing indenta- 
tion load-displacement data is outlined which addresses 
these problems. We begin by presenting data for a variety 
of materials to show that unloading data are usuaUy 
nonlinear. An analysis technique is then presented which 
accounts for the curvature in the unloading data and 
provides a physically justifiable procedure for determin- 
ing the depth that should be used in conjunction with 
the indenter shape function to establish the contact area 
at peak load. The hardnesses and moduli of several 
materials are then computed using the analysis and 
compared with values determined by independent means 
to establish the accuracy of the method, in addition, 
several practical issues concerning procedures for taking 
data and methods for determining load frame compliance 
and indenter shape functions are discussed. 
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il. EXPERIMENTAL 
A. 



The materials used In the study were chosen to span 
a wide range of hardness and modulus. They included a 
soft fee metal (aluminum), a hard bec metal (tungsten), 
two amorphous glasses (soda-lime glass and fused sil- 
ica), and two crystalline ceramics (sapphire and quartz). 
Details concerning purity, finish, and crystal orientation 
arc given in Table L One seemingly obvious choice 
which does not appear in the table is silicon, which was 
intentionally excluded because its indentation behavior 
is complicated by cracking and pressure-induced phase 
transformations, 21 " 0 

Aluminum and tungsten were chosen in preference 
to other metals because both are relatively isotropic 
in their elastic properties; their anisotropy ratios ate 
1.21 and 1.00, respectively? 1 The question of whether 
the modulus measured in an indentation test represents 
that of some specific crystallography direction or some 
average value is therefore not an issue. The same is true 
of the amorphous glasses, which are isotropic in both 
their elastic properties and their hardnesses. Quartz and 
sapphire, on the other hand, have complex hexagonal 
crystal structures, and because of this, results from these 
materials can be used to provide some insight into the 
importance of elastic anisotropy. 

B. Indentation procedure 

All experiments were performed using a Nano- 
indenter* (Nanoindenter is a registered trademark of 
Nano Instruments, Inc., KnoxviHe, IN) at the Oak Ridge 
National Laboratory, a schematic fllustrarion of which 
is Shown in Fig. 2. The system has load and displace- 
ment resolutions of 03 fiH and 0.16 nm, respectively. A 
Berkovich indenter, a three-sided pyramid with an area- 



TABLE I. Summary of 


materials used ra indentation study. 


Material 


Description 


Aluminum 


99.995% pure single crystal mechanically 




polished with colloidal silica 


Tungsten 


99.95% pure single crystal mechanically 




polished with colloidal silica 


Soda-time glass 


Commercial microscope slide 


Fused silica 


Optically 0at substrate material" 


Quartz 


(001) single crystal: optically flat 


Sapphire 


99-995% AJaO*; (001) single crystal 




mechanically polished to optical flatness 




FIG. X A schematic representation of the experimental apparatus 
used to perform the indentation experiments: (A) sample: (B) indenter; 
(Q load application coO; (D) mApiratioit column guide spring 
(E) capaemve A '~ mi — 



to-depth function which is the same as that of a Vfckers 
indenter, was used in all experiments. 

Most experiments were performed using a toad-time 
sequence like that shown in Kg. 3. The indenter was 
first loaded and unloaded three times in succession at 
a constant rate of loading with each of the unloadings 
terminated at 10% of the peak load to assure that 
contact was maintained between the specimen and the 
indenter. The reason for performing multiple loadings 
and unloadings was to examine the reversibility of the 
deformation and thereby make sure that the unloading 
data used for analysis purposes were mostly elastic. 
Affcer die third unloading, the load was held constant 
for a period of 100 s at 10% of the peak value while 
die displacement was carefully monitored to establish 
the rate of displacement produced by thermal expansion 
in the system. Even though the system is thermally 
buffered from its surroundings and the room in which 
IX is housed is temperature-controlled to within *1 °C, 
smalt thermal fluctuations cause some of the machine 




200 300 
TIME (sec) 
FIG. 3. A typical toad-time sequence; peak toad 
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components to expand* ami the expansion is manifested 
in the data as an apparent displacement in die specimen. 
Thermal drift becomes particularly important for small 
indentations made over long periods of time, as was 
the case for some of the indentations in (his study. To 
account for thermal drift, the rate of displacement was 
measured during the last 80 s of the bold period, and 
the displacement data were corrected fay assuming that 
this drift rate was constant throughout die test. In this 
procedure, it is implicitly assumed that the displacements 
measured at the low load during die hold period are 
due mostly to thermal expansion rather than to time 
dependent plasticity in the specimen. Following die bold 
period, the specimen was loaded for a fourth and final 
time, with another 100 s hold period inserted at peak 
load lo allow any final time dependent plastic effects to 
diminish, and the specimen was folly unloaded. 

For each material, six separate peak loads were in- 
vestigated starting at 120 mN and Successively reducing 
the load by a factor of 3 to produce indentations at 120, 
40, 13.3, 4.4, 1*5, and 0-5 mN. The loading/unloading 
rate was also reduced by a factor of 3 starting at a value 
of 5000 pN/s for the 120 mN indentations. The reason 
for reducing the loading rate in proportion to the peak 
load was to assure that all load-displacement curves had 
about the same number of data points since (the system 
samples data at a constant rate of approximately 1.6 
data points per second). Five indentations were made 
at each load, with most of the results presented here 
representing averages for die group. Wherever possible, 
scatter bars representing ± one standard deviation have 
been included in the data. 

The indentation system also had the ability to contin- 
uously measure contact stiffness during indentation This 
is accomplished by superimposing a small oscillation on 
the force signal and measuring the displacement response 
at the same frequency. Continuous stiffness measurement 
was used to provide direct evidence that the contact stiff- 
ness changes continuously during loading and unloading. 
Details of the technique and analysis procedure, parts of 
which have been published elsewhere, 32 '" 34 are included 
in Appendix I. 

111. EXPERIMENTAL RESULTS 

A. General characteristics of load- 
displacement curves 

Wfe begin the discussion of experimental results 
with an overview of the characteristics of the -{pad- 
displacement curves for the six materials studied* Figures 
4-9 present experimental data for each of the materials 
for indentations made to peak loads of 120 mN. The 
differences in hardness of the materials are apparent from 
the large differences in the depth attained a maximum 
load. The softest material is aluminum, with a peak depth 
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FIG. 4. Load versus usdemer displacement for the highest peak toad 
experiment performed on a mechanically polished aluminum single 
crystal. 

of almost 5000 nm, while the hardest is sapphire, which 
was penetrated to a depth of only about 500 nm Both 
the aluminum and tungsten data are typical of materials 
in which the hardness is relatively small compared io 
the modulus, as is observed in most metals; most of 
the indenter displacement in these two materials is 
accommodated plastically, and only a small portion 
is recovered on unloading. The other materials show 
varying degrees of elastic recovery during unloading, the 
largest being that for fused silica. 

The unloading/reloading behavior of the materials 
exhibits some variability. Figure 10 shows an expanded 
view of the sapphire results, demonstrating that the 
unloading and reloading curves tor this material are 
nearly the same. This behavior is also observed in fused 
silica and quartz, with the near perfect reversibility 
suggesting that deformation in these materials after the 
initial loading is almost entirely elastic. 
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FIG. $. toad versus jndemcr displacement for the highest f 
experiment performed <m soda-lime ghaa. 
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200 400 800 *00 1000 1200 
DISPLACEMENT (nm) 
FIG, & Load versus inricntcr tfrspUceroeni for the highest peak load 
experiment performed on fused silica. 



For aluminum, on the other band, the expanded view 
m Fig. 11 shows that the peak load dbplacements shift 
to higher values in successive cycles. In addition, the 
relatively large displacement just prior to final unloading 
is due to creep during the 100 s hold period at peak load. 
Thus, it is apparent that at least some plastk deformation 
occurs in this materia! after the initial loading, and some 
of it is time dependent Hie behavior of tungsten is 
similar to that of aluminum. 

A third type of behavior is thai exhibited by soda- 
lime glass. As shown in Fig. 12, this material exhibits 
distinct hysteresis loops, as might be expected if there 
were a small amount of reverse plasticity upon unload- 
ing. The looping degenerates with cycling, however, and 
after 3 or 4 cycles, the load-displacement behavior is 
largely elastic 

The main reason for discussing the load-displace- 
ment behavior during unloading/reloading is to point out 



that in some materials, displacements recovered during 
first unloading may not be entirely elastic, and because 
Of this, rite use of first unloading curves in the analysis 
of elastic properties can sometimes lead to inaccuracies. 
One way to minimise nonelastic effects is to include 
peak load hold periods in die loading sequence to 
allow time dependent plastic effects to diminish, and 
to use unloading curves obtained after several cycles 
of loading to minimize effects of reverse plasticity. The 
procedure we have adopted here to avoid these problems 
is to load and unload the indentation four tunes and 
include one long period at peak load Our analyses 
are then performed using the final unloading data only. 
Exactly how nonebsttc effects can influence the mea- 
surement of mechanical properties is documented later in 
this paper. 
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FIG. 7. Load versa fcdenter displacement for die highest peak lead experiment performed on a mechanically polished tungsten single 
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the toad versus indenter displacement data for the higb»t peek teed 
oqperitneiit performed on soda-lime glass. 



B. Shapes of unloading curves 

To examine carefully the shapes of unloading curves, 
final unloading data from all six materials are plotted 
On a common set of axes in Fig. 13. For die sake 
of comparison, the displacements have been shifted 
laterally by subtracting the final depth, h fm from the 
total depth, ft, to force all the curves to pass through a 
common origin. Careful examination of the figure shows 
that none of the data is linear; rather, each unloading 
curve is slightly concave up over its entire span. This can 
be better appreciated by plotting the data logarithmically, 
as has been done in Fig. 14. The fact that the data are 
linear on logarithmic axes implies that the unloading 
curves arc well described by a power law relation like 
that of Eq, (1). The exact values of the experimentally 




4300 4400 4500 4600 4700 
DISPLACEMENT (rrm) 



4000 



FIG. 11. An expanded view of the unloading/reloading portion of 
the load vervtu indenter displacement data for the tugW peak k»d 
experiment performed an a rncehanicaiiy polished aluminum stogie 
crystal. 



observed power law exponents, as well as the correlation 
coefficients for the power-law fits, are given in Table II. 
The power law exponents for the six materials studied 
here vary from 1.25 to 131, and the fact that they are 
distinctly greater than one implies that none of the data 
is consistent with flat punch behavior. 

C. Continuous stiffness measurement 

Additional evidence mat the flat punch approxi- 
mation is inappropriate is provided by the continuous 
measurement of stiffness. Here, we discuss results for a 
material not included in Table I, ekctxopolished tung- 
sten. This material is particularly interesting because 
well -annealed single crystals indented by a Berkovich 
indenter at sufficiently tow loads can be deformed purely 
elasticaHy— that is, the contact impression fully recovers 




100 200 300 400 300 GOO 
DISPUU^MENT,h-h # <nm)' 

ITG. 13. The final unloading segment of the load wsus indenter 
displacement data from the highest peak load experiment* for all six 
materials plotted so that the curves pass through a common origin. 
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upon unloading, no hardness impression remains after 
indentation, and the loading and unloading corves trace 
one another perfectly, 3133 As shown in Fig. 15, this 
behavior was observed in electropolished tungsten for 
peak loads of 0.5 mN and less. For purely elastic contact, 
both the co ntact area and the contact stiffness are unique, 
reversible functions of the load, and they increase and 
decrease with load according to simple functional rela- 
tionships that depend on the geometry of the indenter. 
Because of this, contact stiffness correlates directly with 
contact area, and hence, changes in contact stiffness in an 
indentation experiment can be correlated with changes in 
contact area. Continuous stiffness measurement results 
for the indentation data io Eg. 15 are presented in 
Fig. 16, plotted as stiffness versus time. Comparison of 
these data with the load-tune history shown in Fig. 3 
shows that the measured contact stiffness, and thus the 
contact area, does indeed increase and decrease in the 
way char would be expected based on the loading history 
(note that the hold periods in Rg. 3 were not included 
in acquiring the data of Figs. 15 and 16). 

At higher peak loads, the indentation contact to 
electropolished tungsten is not just clastic. Rather, as 
shown in Fig. 17, when a threshold load of about 1.0 mN 



TABLE II- Pur^cters describing power kw fitj of unloading ouves. 



Material 


A 

(mN/wnT) 


m 


Correlation 
coefficient 


Aluminum 


0265 


13* 


0.99993$ 


Qmnz 


OJ0215 


L43 


0.999985 


SodVlimc glass 


0JO279 


1.37 


0599997 


Sapphire 


0.0435 


1.47 


0.999998 


Fused silica 


OjOSOO 


1.25 


0.999997 


Tungsten 


0.141 


«1 


0999986 



0.6 
0,5 

I" 

" 0JJ 
02 

0.1 



■ ' I ' — ' ' 1 1 1 J ' I 

ELASTIC CONTACT 
TUNGSTEN 



20 



25 



0 5 10 15 

DISPLACEMENT (nm) 

HO- 15. toad versus indentcr displacement lot a fully elastic contact 
on an etruiepuldsbed angle crystal of tungsten. The data consist of 
tour sequential taadmr/Junlcadmg cycles. 



is reached, a sudden jump in displacement corresponding 
to the onset of plasticity Is observed, and a permanent 
hardness impression is formed. The continuous stiffness 
measurement for this indentation is shown in Fig. 18. it 
is seen that for each of the four unloading*, the contact 
stiffness changes immediately and continuously as the 
specimen is unloaded. Thus, the contact area, which 
varies in the same way as the contact stiffness, is not 
constant during the unloading of the plastic hardness im- 
pression, even during the initial stages of unloading, and 
the flat punch approximation is therefore questionable. 
Continuous and immediate changes to contact stiffness 
during unloading were also observed in each of the other 
materials examined in this study. Tungsten was chosen 
to illustrate the behavior because of its unique elastic 
behavior at tow loads. 
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FIG. IS. Contact stiffness versus lime for a fully ctaslic contact an an 
ekarupoltthed single crystal of tungsten measured with the continu- 
ous Stiffncsx measurement technique. 
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DISPLACEMENT (ran) 

FIG. 17. Lewd versus itdenter displacement far a plastic contact on 
an electrapoiisbed single crystal of tungsten. The dam consist of four 
sequential loadingAinkwfin^ cycles- 

D. Comments on linear unloading 

Before leaving this section, a comment is warranted 
as to why the observations here are different from other 
reports in the literature which suggest that unloading 
curves in some materials are linear, at least in the initial 
stages of unloading. The primary observations of this 
sort are those of Doemer and Nix, who report linear 
unloading in metals over most of the unloading range 
and in silicon for at least the first one-third of the 
unloading curve. 3 The data they present as representative 
of metallic behavior arc that of aluminum. The data 
consist of one loading and unloading plotted with the 
unloading and loading curves on a single set of axes. On 
such a plot, the unloading curves are so steep that they 
give the appearance of being linear even though they 
may not be. In Fig. 4, for example, the unloading curves 
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FIG. IB* Contact snf&tcaa versus tunc foe a plastic contact on an dec* 
trOfjclehed single crystal of tungsten measured using the 
ftiflness meafiufentenl technique. 



for aluminum appear to be linear, but when plotted on 
expanded or logarithmic axes such as those in Figs. U, 
13, and 14, the nonlinearity becomes apparent. We also 
wish to note that improper values for the load frame 
stiffn e ss can influence unloading curves in a way that 
makes them appear to be more linear than they really are. 

With respect to die observations of linear unloading 
to silicon, we note that the indentation behavior of this 
material stands ina class Of its Own. Studies have shown 
that a reversible., pressure induced transformation to a 
denser phase occurs beneath the mdenter**" 30 When the 
indenter is withdrawn, the transformation reverts, and a 
portion of the unloading displacement is recovered by a 
nonelastic process. Unloading curves for silicon are thus 
unlike those tor any other material, and they sometimes 
may appear to have a linear region near the top. In our 
experience, however, careful examination of such data 
has shown that it is curved wirh power law exponents in 
the range reported here for the other materials. 

IV. A NEW METHOD OF ANALYSIS 

The observations presented in the previous section 
suggest that the flat punch approximation is not an 
entirely adequate description of the behavior of mate- 
rials when indented by a Berkovich indenter. Here, we 
propose a new method of analysis based on analytical 
solutions for other indenter geometries. In addition to 
accounting for curvature in the unloading data, the 
method also provides a physically justifiable procedure 
for determining the depth that should be used in con- 
junction with the indenter shape function to establish 
the contact area at peak load. 

As mentioned in the introduction, Sneddon has de- 
rived closed form analytical solutions for punches of 
several geometries. In addition to the flat punch, he 
has also considered punches with conical geometries 
and paraboloids of revolution.^ Wc proceed here on 
die assumption that the behavior of one of these lat- 
ter geometries gives a better description of the elastic 
unloading of an indentation made with a Berkovich 
indenter. 

The conical indenter is a natural choice since, like 
the Berkovich indenter, its cross-sectional area varies as 
die square of the depth of contact and its geometry is 
singular at the tip. The paraboloid of revolution, whose 
behavior is the same as that of a spherical indenter in 
the limit of small displacements, is also a potentially 
useful geometry in that no real indenter is ever perfectly 
sharp; i.e., at some scale the tip exhibits some rounding. 
In addition, because of plasticity, an clastic Singularity 
cannot really exist at the tip of the indenter, and the 
paraboloid geometry potentially accounts for this. For 
both geometries, the load-displacement relationships are 
nonlinear and the contact area changes continuously 
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during unloading. An analysis of the geometry that 
works best is presented in the next section by compar- 
ing predictions of the method with actual experimental 



An important fundamental assumption in the de- 
velopment is that the equations describing the elastic 
un loading of a fiat, serin-infinite half space are the same 
as those for an indented surface; that is, Sneddon's 
solutions apply equally well Id a fiat surface or a surface 
with a hardness impression. The justification for this as- 
sumption is based largely on die empirical observations 
discussed in the introduction that hardness impressions 
formed in metals by conical and spherical tndentcrs are 
also conical and spherical with slightly different indnded 
op angles and radii As discussed in greater detail 
elsewhere, 25 Che mathematical description of the elastic 
loading and unloading of these perturbed surfaces is 
exactly the same as that for the fiat surface if adjustments 
are made to the geometrical parameters describing the tip 
angle of die cone and the effective radius of the sphere. 
Since spherical and conical indentera represent two very 
different geometries, Le_, one very sharp and one very 
blunt, it seems reasonable that this behavior may hold 
for other axisymmetric indenters as well. 

Figure 19 shows a cross section of an indentation 
and identifies the parameters used in the analysts. At 
any rime during loading, die total displacement h is writ- 
ten as 



<5) 



where h c is the vertical distance along which contact 
is made (hereafter called die contact depth) and fa is 
the displacement of the surface at the perimeter of the 
contact. At peak load, the load and displacement are 
Ptntx and Aaax, respectively, and the radius of the contact 
Circle is a. Upon unloading, the elastic displacements are 
recovered, and when the indenter is fully withdrawn, the 
final depth of die residual hardness impression is */. 




The experimental parameters needed to determine 
hardness and modulus are shown in the schematic load- 
displacement data shown in Fig. 20. The three key 
parameters are the peak load die depth at peak 
load (htttaX and the initial unloading contact stiffness 
It should be noted that the contact stiffness is 
measured only at peak hnd» and no restrictions are 
placed on the unloading data being linear during any 
portion of the unloading- 

The analysis begins by rewriting Eq (3) as 



2 yfA 



(6) 



which relates the reduced modulus, to the contact 
area, A, and the measured stiffness, 5. As discussed 
previously, this relationship holds for any indenter that 
can be described as a body of revolution of a smooth 
function and is thus not limited to a specific geometry. 25 
Measurement of the initial unloading slope can thus be 
used to determine the reduced modulus if the contact 
area at peak load can be measured independently. 

The area of contact at peak load is determined by 
the geometry of the indenter and the depth of contact, h c . 
Following Oliver ex cL, ,a we assume thai the indenter 
geometry can be described by an area function F(h) 
which relates die cross-sectional area of the indenter to 
the distance from its tip, k Given that the indenter does 
not itself deform significantly, die projected contact area 
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FIG. 19. A schematic itpreseaUlkM of a section through 
Uqq showing various quantities USCd in the analysis. 
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FIG. 20- a. schematic repiesentanon of kwd versus indenter displace- 
ment showing qcuntmes toed is the analysis as weu as a graphical 
tmerpfetatioa of the contact depth. 
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at peak load can then be computed from the relation 

A«F(ft c ) (7) 

The functional form of F must be established expert- 
mentally prior to analysis. 

To determine the contact depth from die experimen- 
tal data, it is noted that 

h< = Knsx — (8) 

which follows directly from Eq, (5). Since can 
be experimentally measured, the key to the analysis 
then becomes how the displacement of the surface at 
the contact perimeter, can be ascertained from die 
load-displacement data* 

The deflection of the surface at the contact perimeter 
depends on die indenter geometry. For a conical indenter, 
Sneddon's expression for the shape of die surface outside 
the area of contact 16 can be used to give 



The quantity (A — hf) appears in this expression rather 
than h by itself since Sneddon's solution applies only 
to the elastic component of the displacement. In addi- 
tion, Sneddon's face-displacement relationship for die 
conical indenter yields 

(h-h f )^2~ (10) 

where 5 is the stiffness. Substituting Eq. (10) into Eq. (9) 
and noting that the contact area of interest is that at peak 
load, one obtains 

A. = <%^ (11) 

where the geometric constant £ for the conical indenter 
is given by 



(12) 



or € = 0-72. When similar arguments are made for the 
flat punch and the paraboloid of revolution, Eq. (11) 
is again obtained with different geometric constants. 
For the flat punch, e — 1. and for the paraboloid of 
revolution, e = 0.75. 

The graphical interpretation of Eq. (11) is shown 
in Fig. 20. For € - 1, the value for the flat punch, 
h, = Pjw/S, and the contact depth h e is given by 
the intercept of the initial unloading dope with the 
displacement axis, interestingly, this is precisely the 
depth used by Doeroer and Nix in their analysis based on 
the fiat punch approximation. Urns, the current method 
is consistent with the Doerner and Nix approach when 
the fiat punch geometry is assumed. For the conical and 
paraboloid indenters, however, the contact depths art 



greater than those for the flat punch, and this must he 
accounted for in analyses using these indenter geome- 
tries if accurate measurements are to be obtained. The 
range of *r for the tndenters considered here is shown 
in Fig. 20. 

In adoption to the modulus, the data obtained using 
the current method can be used to determine the hard- 
ness, H We define the hardness as the mean pressure the 
material will support under load. With this definition, the 
hardness is oomputcd from 



f*\ max 



(13) 



where A is die projected area of contact at peak load 
evaluated from Eq. (7). It should be noted that hardness 
measured using this definition may be different from 
that obtained from the more conventional definition in 
which die area is determined by direct measurement of 
the size of the residual hardness impression, The reason 
for the difference is that, in some materials, a portion 
of the contact area under load may not be plastically 
deformed, and as a result, the contact area measured 
by observation of the residual hardness impression may 
be less than that at peak load. However, for most of 
the materials considered in this study, the contact areas 
computed using the procedure outlined here compare 
favorably with residual contact areas measured in the 
SEM, and the two definitions of hardness give similar 
results. This will be discussed further in the next section. 

V. IMPLEMENTATION OF THE METHOD 
A. Choice of indenter geometries 

The first step in implementing the new analysis 
procedure is to identify which of the indenter geometries 
best describes the experimental data and therefore which 
value of t should be used While it is conceivable that die 
best description varies from one material to another, the 
data presented in Fig. 21 suggest otherwise. The figure 
shows the unloading curves of Fig. 13 replotted with the 
loads and displacements normalized with respect to their 
maximum values. When plotted this way, it is apparent 
that the unloading data are all remarkably similar in 
form In fact, with the exception of fused silica data, 
which are slightly higher than the rest, the curves for 
the materials are so similar that they can hardly be 
distinguished. This suggests that unloading behavior of 
a wide variety of materials can be described by a single 
indenter geometry. 

As for which indenter geometry works the best, 
die primary due is found in the power taw exponents, 
nu describing the unloading behavior. Table 111 summa- 
rizes the power law exponents predicted by Sneddon's 
analyses along with the associated values of c The 
experimentally determined values of m listed in Table 11 
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FIG. 21. Normalized final unloading carves tor ail six 
showing similarities in ahape. 



range from 1.25 to 1-51 with a mean value of 1.40. 
Comparison of Che two tables thus shows thai based 
on the power law exponents, the unloading behavior 
is best described by the paraboloid geometry, and the 
value of € we use in the analysis procedure is 0.75. Why 
the paraboloid is a better description of the unloading 
behavior is probably due to the fact that the elastic 
singularity characteristic of the conical geometry is not 
physically realizable when plasticity occurs, and the 
pressure distribution which actually forms around the tip 
of the indenter is more lite that predicted by the parabola 
of revolution. It is notable, however, that even had the 
conical geometry been chosen, the only difference In the 
analysis would have been to use a slightly smaller value 
of c, ie., € - 0.72. 

B. Measurement of the Initial unloading stiffness 

One important practical question that arises in the 
analysis procedure is how the initial unloading stiffness, 
S, should be measured from the unloading data. One 
simple way to accomplish the measurement is to fit a 
straight line to a fraction of the upper portion of me 
unloading curve and use its slope as a measure of the 
stiffness. The problem with mis is that for nonlinear 
unloading data, the measured stiffness depends on how 
much of the data is used in the fit 

To illustrate how important this can be, plotted 
in Fig. 22 arc stiffnesses measured for tungsten as a 

TABLE HI. Punch parameters used in data analyses, _ 



Punch geometry 
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0.7S 
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FfG. 22. The peak toad stiffness of a 120 mN indentation in tungsten 
computed using linear fits of the unloading data and the power law 
fitting technique. The dashed tines are power law fits applied to the 
first and last unloading segments. 



function of the fraction of the unloading data used 
in the fit Two features are worthy of note. First, the 
stiffnesses obtained from the first and last unloading* 
are very different. The reason for this is that there is a 
significant amount of creep during the first unloading, 
causing the slope of the upper portion of the unloading 
curve to be abnormally high. Because of this, mere is no 
way that the proper unloading stiffness can be identified 
from the first unloading data using this technique. These 
effects can be minimized by the inclusion of peak load 
hold periods in the loading sequence to arminish time 
dependent plastic effects, as was done just prior to the 
final unloading in this study. Second, it is seen mat even 
in the final unloading data, there is a significant variation 
in the measured stiffness, depending on how much of 
the unloading curve is used in me fit. For this reason, an 
alternative procedure is desirable. 

The procedure wc have adopted is based on the 
observation that unloading data are well described by a 
simple power law relation, as directly evidenced by the 
excellent correlation coefficients in Table II. The actual 
relationship we use to describe the unloading data for 
stiffness measurement is 



p = A{h-h,y 



(14) 



where the constants A, m, and hj are all determined by 
a least squares fitting procedure. Hie initial unloading 
slope is then found by analytically differentiating this 
expression and evaluating the derivative at the peak load 
and displacement 

Results for tungsten are included in Fig. 22 as 
horizontal dashed fines. It is seen that the stiffness 
computed from the final unloading curve is very close 
to thai of the linear firs obtained using small fractions 
of the unloading data. Interestingly, when applied to the 
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first unloading data, the power law method produces a 
stiffness that is only slightly greater ten that derived 
from the last unloading curve, thereby implying that the 
power law method is kss sensitive to creep. For these 
reasons, the power law method is preferred and is used 
in all subsequent analyses. 

C« Determination of toad frame compliance 
and diamond area function 

Another important practical consideration is how 
well the load frame camp] iaitce and area function of the 
diamond are known. The load frame compliance is im- 
portant because the measured displacements are the sum 
of the displacements in the specimen and the bad frame, 
and, thus, to accurately determine specimen displace- 
ments, the load frame compliance must be known with 
some precision. This is especially important for Urge 
indentations made in materials with high modulus for 
which the load frame displacement can be a significant 
fraction of the torn! displacement. The area function is 
equally important since computation of both the modulus 
and the hardness depend on the contact areas through 
Eqs. (6) and (13). 

The standard procedure used in the past tor de- 
termining area functions has been to make a series of 
indentations at various depths in materials in which 
the indenter displacement is predominantly plastic and 
measure the size of the indentations by direct imaging. 
Using the contact depths computed from the indenta- 
tion load-displacement data, an area function is then 
derived by empirically fitting a function to a plot of 
the imaged areas versus the contact depths. Clearly, the 
imaging technique must have sufficient contrast, reso- 
lution, and precision to accurately determine the areas. 
Optical imaging works well for larger indentations, but 
cannot be used exclusively since it is usually necessary 
to characterize the area function well into die submicron 
range. Scanning electron microscopy is also of limited 
value for small indentations because the shallowness of 
the hardness impressions results in very poor contrast 
TEM replication methods have proved useful, but they 
are tedious, time-consuming, and require equipment and 
expertise which may not be available in all laboratories. 

Here, we propose a method for determining area 
functions that requires no imaging at all. The method 
is based on the assumption that the elastic modulus 
is independent of indentation depth, justification for 
which will be given shortly when results of contact 
area measurements for several different materials are 
presented. It is important to note that specific values for 
the modulus arc not assumed. 

The method follows by modeling the load feme and 
the specimen as two springs in series, in which case 



where C is the total measured compliance, C, is the 
compliance of the specimen, and Cf is die compliance of 
the load frame. Since the specimen compliance during 
elastic contact is given by the inverse of the contact 
stiffness, 5, Eqs. (3) and (15) combine to yield 



(16) 



It is thus seen that if the modulus h constant, a plot of C 
vs A* 1 is linear for a given material, and the intercept of 
the plot is a direct measure of the load frame compliance. 
The best values of C f axe obtained when the second term 
on the right-hand side of Eq. (16) is small, Lc., tot large 
indentations. 

To find the area function and the load frame com- 
pliance, wc take advantage of the fact that relatively 
large indentations can be made in aluminum because of 
its low hardness. In addition, for the largest aluminum 
indentations, the area function for a perfect Berkovich 
mdeoter 



(17) 



c = c a + c f 



(15) 



can be used to provide a first estimate of the contact 
area. Initial estimates of Cf and E r were thus obtained 
by plotting C vs A Aa for the two largest indentations 
in aluminum. Using these values, contact areas were 
then computed for all six indentation sizes by rewriting 
Eq. (16) as 

a * 1 1 

from which an initial guess at the area function was made 
by fitting the A vs data to the relationship 

A(h c ) _ 24.5n* + C 2 fti + C 2 hl' 2 + Ca^ /4 

+ + C*^ /U * (19) 

where Ci through C 8 are constants. The lead term de- 
scribes a perfect Berkovich indemcr; the others describe 
deviations from the Berkovich geometry due to blunting 
at the tip. 

The procedure is nor complete at this stage because 
the exact form of the area function influences the values 
of C f and E rr so using the new area function the 
procedure was applied again and iterated several times 
until convergence was achieved. Figure 23 shows a plot 
of die final values of (C - C f ) vs A*** m Note mat the 
data are linear and extrapolate to (C - C/) = 0, as they 
should once die proper load frame compliance and area 
function have been found. 

To extend the area function to smaller depths and to 
check the validity of the constant modulus 'assumption, 
the procedure was subsequently applied to all the other 
materials, in doing so, the load frame compliance was 
held constant at the value determined from the aluminum 
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FIG. 23. A plot of <C - Q) vsA* 4 ^ for atommwn. Hid error btn are 
two standard deviations in length. 

data since the largest contacts were made in this ma- 
terial. The area of Hie largest indentation in each of 
the other materials was computed using die aluminum 
area function* and die modulus for each material was 
estimated using Eq. (6). Hie areas needed to produce a 
constant modulus as a function of depth for each material 
were then computed using Eq. (18> The resulting areas 
and the corresponding contact depths for all materials 
were fit as a group using Eq. (19) to define a new 
area function. The process of detennining E r at the 
largest depth for each material (using the latest area 
function), recalculating areas assuming contant modulus, 
and redetermining the area function was then repeated 
iterativeiy until the area function accurately described 
the data for all of die materials at all depths. Note 
that nowhere in the procedure was a specific value for 
a modulus assumed; the only assumption is that the 
modulus is independent of depth. 

The results for all the materials arc plotted in Fig. 24 
as contact area versus contact depth. The solid line 
through the data is the final composite area function 
derived by fitting Eq. (19) to all die data. The fact that 
one single area function adequately describes the data 
suggests that the constant modulus assumption and the 
use of a single area function are valid. The composite 
area function derived from these data is used in all 
further analyses. 

Subsequent to this work, a more specific procedure 
was developed for detennining area functions and 
load frame compliances using aluminum and fused 
quartz as calibration materials. Details are described in 
Appendix fl. " s 

D. Assessment of predictive capability 

To assess the predictive capability of the method, 
load-displacement data for each of the six experimental 
materials were analyzed according to the procedures 
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described in Sec IV so determine contact areas, moduli, 
and hardnesses. The elastic constants of the diamond 
indenter used m Eq- (2) to compute the specimen modu- 
lus from the reduced modulus were £ f = 1141 GPa and 
i/i = 0.07, Here, we compare predictions of the method 
with values assessed using independent means. 

The first comparison is between the contact areas 
calculated from the load displacement data and contact 
areas measured by direct imaging in the SEM. Figure 25 
presents results for all six materials plotted as calculated 
versus imaged areas. The good agreement between the 
two measurements suggests that the method works wcl! 7 
with die only possible exception being aluminum at 
high loads, for which the calculated areas arc slightly 
greater than imaged areas. That aluminum behaves in 
this way may be related to the fact that it was the only 
material Studied here for which there was significant 
pile-up around the indentation. 

The good agreement in Fig. 25 also suggests that the 
peak toad contact areas computed using the new method 
and the imaged residual contact areas are much the same 
in most materials, and thus the definition of hardness 
given in Eq. (13) is consistent with the conventional 
definition, litis finding was at first somewhat surprising 
to us, since there is significant elastic recovery upon 
unloading in the ceramics and glasses we examined (see 
Figs. 5-8)> and we expected to find the residual contact 
areas in these materials to be measurably smaller than 
the peak load values. The reason for the agreement can 
be at least partially understood by examination of the 
SEM micrograph In Fig. 26. The unusual shape of this 
fused silica indentation suggests that during unloading, 
the Sides of the indentation elasticaily recovered white 



J. Motor. RfiA, Vol. 7, No. 6, Jurt 1092 



1S77 



TThrrcQOTrrc n i 



JUN 22 2005 14:02 FR CISTI ICiST 



TO 18027696501 



P. 17 



W, C. Offver and G> M_ Pharn An improved technique for tfatermimng hardness ami elastic modulus 



10 3 



* 6 1<r* 

a. 



1Q 1 



0.1 



• • 






CONTACT AREA 






: CORRELATION 








o 


ALUMINUM 




o 


QUARTZ 




o 


SODA UMS GLASS 




X 


SAPPHIRE ~ 












TUNGSTEN 


1 I 







0.1 



1000 



RC 



1 10 100 

WAGED AB€A (lim 1 ) 
25. Contact areas computed using the new analysis procedure 
nUct areas for ail six materials. The enor tar? an 
Standard deviations in length. 



the corners did not, thus producing the odd appearance 
of the edges Of the impression. Apparently, the additional 




FIG. 26. A scanning electron micrograph of a 4G mN 
fused silica. 



plasticity caused by the stress concentration at the edges 
the indemer serves to mark permanently the position 
of the corners of the indentation at peak load Thus, 
if the area of the hardness impression b measured 
by computing the area of the triangle defined by the 
imp l&ssnm comers, as is most often done in convt^itional 
hardness testing and as was done here, then the size of 
the imaged indentation is a good measure of the peak 
toad contact area, la thus regard, it is notable that the 
same agreement would not be expected for indentations 
made with smooth, axiatly symmetric tndeoters (Le. t 
cones, paraboloids or spheres), for which elastic recovery 
would occur uniformly around the periphery of the 
impression. 

The hardnesses computed from the calculated areas 
are shown in Fig. 27. The data show that there is very 
little indentation size effect (ISE) in any of the materials 
with the exception of aluminum and tungsten in which 
there is a modest increase in hardness at low loads. This 
could be due to surface-localized cold-work resulting 
from polishing. 

The moduli predicted from the indentation load- 
displacement data are shown in Fig. 28. Note thai to 
compute the modulus from the measured values of E r 
requires chat Poisson's ratio be known. The assumed val- 
ues for Potsson's ratio art shown in Table IV, He data in 
Fig. 28 show that once again, there is very little evidence 
for an indentation size effect; U., the moduli remain 
more or less constant over the entire range of load. 

To assess the predictive capabilities^ the modulus 
determination, moduli obtained by averaging the experi- 
mental results at the two highest loads are compared 
with commonly accepted values from the literature in 
Table IV. For the crystalline materials, the literature 
values were obtained by averaging the Vbigt and Reus* 
average elastic constants (or the Hashin and Shtrikman 
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average elastic constants for the cubic materials) using 
the data tabulated 1 in Simmons and Wang.* 5 

Hie table shows that for aluminum, tungsten, soda 
lime glass, and fused quartz, the moduli computed from 
the indentation load-displacement data are all within 4% 
of the literature values. Since these materials ail exhibit 
very isotropic elastic properties, there is no question as 
to what the measured modulus should be, and the good 
agreement between the experiment and literature moduli 
demonstrates that the method proposed here works welL 
For sapphire and quartz, on the other hand, the agreement 
in the table is not quite as good. Hie measured modulus 
is higher than the Mngt/Reuss average by about 9% for 
the sapphire and about 30% for the quartz. However, 
both these materials ait highly anisotropic, and because 
of this tt t$ not clear to what value the measured modulus 
should be compared. Ou the one hand, h can be argued 
that the most appropriate modulus is that in the direction 
of testing, since the elastic displacements are primarily 
in this direction. On the other hand, since the formation 
of the contact impression involves deformation in many 
directions, it can be argued that the measured modulus 
should be some average quantity. 



Unfortunately, which of these hypotheses is more 
likely, if either, cannot be determined using the limited 
data obtained in this study. For sapphire, it could be 
argued that the reason that the Votgr/Renss average of 
403 GPa is tow is because the measured modulus of 
441 GPa is strongly influenced by die 499 QPa modulus 
in the direction of testing (the oauris). Hie same, 
however, cannot be argued for quartz since the experi- 
mentally measured modulus (124 GPa) is greater than 
either the Vbigt/Rsuss average (95.0 GPa) or the c-axis 
modulus (105 GPa)* The influence of elastic amsotropy 
on the measurement of modulus using indentation 
load-displacement methods is an area that requires 
further study. 

VL CONCLUSIONS 

(1) Careful examination of Indentation load- 
displacement data obtained for six materials using a 
Berkovicfa indenrer reveals that unloading curves can be 
accurately described by the power law relation 

P=A(h-h f ) m 

where P is the load, (h-h f ) h the elastic displacement, 
and A and m are material constants. The fact that the 
power law exponents in the above relation are always 
greater than 1 (they vary from material to material in the 
range 1-25 to 1.51) suggests that the flat punch method 
of analysis for determining hardness and modulus from 
indentation load-displacement data is not entirely ade- 
quate. This is corroborated by dynamic measurements of 
contact stiffness which show that contact area changes 
continuously as the indenrer is withdrawn. 

(2) The curvature in the unloading data can be 
accounted for by assuming that the mdenter behave* 
as a punch whose geometry is such that the contact 
area changes continuously during unloading. Conical 
and paraboloid punches satisfy this condition and arc 
attractive for two other reasons: (1) analytical solutions 
for the elastic contact problem are available for each, 
and (2) each has features which model die Berkovich 
geometry in different ways. Careful examination of the 



TABLE IV. Comparison of measured moduli with values in the literature 





Experimental 


Standard 


Litetature 
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modulus (GPa) 
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shape of unloading corves suggests that the paraboloid 
is the preferred geometry. 

(3) Using a new analytical procedure, contact areas, 
hardnesses, and moduli have been computed from the 
indentation load-displacement data for a variety of peak 
loads. Comparison of predicted contact areas with values 
measured independently in the SEM shows that the 
contact area determined by the analysts procedure is 
a good measure of the size of Ac residual contact 
impression. As a consequence; there is good correlation 
between hardness measured from the load-displacement 
data and that measured using conventional techniques* 
The moduli predicted by the method are within 4% of 
values reported in the literature for the materials with 
isotropic elastic properties. For anisotropic materials, a 
question stilt remains as to what die modulus measured 
by load and displacement sensing indentation techniques 
represents. 

(4) The procedure works well only when several 
key quantities are known with some accuracy and pre- 
cision. These include die unloading stiffness, the load 
frame compliance, and the area function of the diamond. 
Procedures for determining each are outlined. The pro- 
cedure for determining me area function is attractive 
because it docs not require imaging of Indentations, 
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APPENDIX I: Continuous measurement of 
contact stiffness by a dynamic technique 

The continuous measurement of stiffness during 
indentation is accomplished by applying a small oscil- 
lation to the force signal at a relatively high frequency 
(693 Hz). The amplitude of the force oscillation is kept 
sufficiently small that the deformation process is act 
significantly affected by its addition. The corresponding 
displacement oscillation is monitored at the excitation 
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frequency using a two-chanael, phase-sensitive detec- 
tor (lock-in amplifier) which can accurately measure 
displacements as small as 0.001 nm using frequency 
Specific amplification. The phase difference between the 
displacement signal and the force signal can also be 
measured A time constant of 0.33 s provides a good 
combination of low noise and dynamic response lor the 
results described in this work. 

Determination of the stiffness of the contact requires 
an understanding of the response of the entire system, 
a dynamic model for which is shown in Fig- AL An 
analysis of this model reveals that the stiffness of the 
contact 3, can be calculated from the amplitude of the 
displacement signal from 

(Al) 

or from the phase difference between the force and 
displacement signals from 

tan(4)r= (A2> 

where 

C f = the compliance of the load frame (^-1.13 m/MN) 
Kj =* the stiffness of the column support springs 

<~«0 N/m) 
D = tiie damping coefficient (~S4 N s/ro) 
P„, = die magnitude of the force oscillation 



DYNAMIC MODEL 




Mass^m 



FIG. Al. A dynamic nuxicJ for tfae system. 

J. Mater. Re*.. 



h(J) — die magnitude of die resulting displacement 
oscillation 

w = frequency of the oscillation (693 Hz) 

= the phase angle between the force and displace- 
ment signals 
m — mass (—4-7 gms) 

The constants that must be known fee die analysis 
are Q, and D. The constant 1^ and D are experimen- 
tally determined by examining the motion of the system 
according to die model in Fig. Al when the tndenter 
is not touching the surface (5 = G)_ The technique for 
determining K m is an important part of the analysts 
technique described in the body of this work and is 
fully described there. Details of the damping coefficient 
measurement and the calibration of die system have been 
described elsewhere. 33 - 34 

For a complete understanding of the technique, 
several characteristics of the system should be noted. 
First, due to the double gap in the displacement sensing 
capactrive gage, the damping coefficient D varies as the 
fourth power of the displacement of the moving plate and 
is thus not constant for large plate motions. However, 
since the positional derivative of D near the center of 
travel is small, D is relatively constant over the range of 
positions used in a normal indentation experiment. 

Second, by defining a reduced spring constant 
K = (S" ! + CfY x -h iC, the model can be reduced to the 
simpler system shown in Fig- A2. The natural frequency 
for this system is 




(A3) 



Since the critical damping coefficient is 4 Nm/s. the 
system Is over-critkally damped. This means that the 
displacement amplitude is single valued and no resonant 
peak is observed. Typically, the characteristics of the 
components of an oscillating system ate fixed and the 
system response is examined by changing the frequency 
of oscillation. In die case considered here, the excitation 
frequency is held constant while the stiffness of the 
reduced spring constant K is increased from a very small 
value = 50 N/m) to a high value (5 x 10? N/m), 
This means initially the system is operating above the 
resonant frequency and the phase angle is greater than 
90 degrees. As the stiffness S, and therefore K, increases, 
the resonant frequency of the system increases and ap- 
proaches the excitation frequency. At the highest values 
of K measured in this work, the resonant frequency of 
the system is greater than the excitation frequency and 
the phase angle approaches zero. As has already been 
noted, even though the system goes through resonance, 
there is no peak in die displacement amplitude due Co the 
high damping coefficient. Clearly the complete solutions 
for the motion of the system given in Eqs. (Al) and 
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Mass = rn 



FK3. A2. A simpler model for das system whictl 
the spring 



oy redefining 



(A2) must be used with no simplifying assumptions. 
In practice, the value of 5 calculated from the phase 
signal works well only at small stiffnesses where die 
phase angle is of some magnitude At high stiffnesses 
the phase angle becomes less than a few degrees and is 
too small relative to the noise. The value of stiffness 
calculated from the amplitude measurement is more 
accurate, particularly at higher stiffnesses. 

APPENDIX B. indent* ar*a function 
calibration procedure 

The area function of die tndenter relates its cross- 
sectional area to die distance from its tip. In the course 
of performing this study, a procedure was developed for 
determining the area function which docs not require 
the imaging of indentations, The procedure involves 
making a series of indentations in two standard mate- 
rials — aluminum and fused quartz — and relies on die 
facts that both these materials are elasticaily isotropic, 
their moduli are well known, and their moduli arc 
independent of indentation depth. 

The procedure requires that indentation contact stiff- 
nesses be measured accurately and precisely over a- wide 
range of indentation depths. The first step in determining 
these stiffnesses is to determine precisely the load frame 
compliance. This is best accomplished by indenting a 
material in which large, very stiff contacts can be made at 
relatively low toads. We have found dial wcti-anncaled, 
high purity aluminum is an excellent choice: it is readily 



available, has a low hardness, and is nearly etasttcally 
isotropic. 

Some care must be exercised in preparing the alu- 
minum to assure that its sorraee is Smooth and unaffected 
by work hardening. The mechanical polishing procedure 
we have adopted is as follows: (!) Mount the specimen in 
epoxy resin to form a standard metallurgical mount with 
a small section of heavy-walled stainless steel tubing 
surrounding the aluminum to prevent plastic straining 
when clamping forces are applied to the mount during 
testing. (2) Grind the specimen with Successively finer 
grits of SiC abrasive paper through 600 grit, being 
sure to remove the damage layer at each step. (3) 
Polish for 8- 12 h on a vibratory polisher using 3 $aa 
diamond paste and water on a TEXMET™ doth fol- 
lowed by 8-12 h using 03 jrm diamond and water on a 
MASTERTEX™ doth. (4) Final polishing is performed 
using colloidal silica and a MASTERTEX™ cloth for 
15-20 min. 

A series of indentations are made in the aluminum 
using the fust six peak loads and loading rates shown 
in Table AL The load time history we recommend is 
as follows: (1) approach and contact surface, (2) load 
to peak load, (3) unload to 90% of peak load and hold 
for 100 s, (4) reload to peak load and hold for 10 s, 
and (5) unload completely at half the rate shown in 
Table AL The lower hold is used to establish thermal 
drift and the upper held to minimize time-dependent 
plastic effects. The final unloading data are used to 
determine the unloading stiffnesses using the power law 
fitting procedure described in Sec. V.B. 

The load frame stiffness is determined from the 
aluminum data using an iterative procedure similar to 
that outlined in Sec. V.C The procedure is actually 
somewhat simpler because the known elastic constants 
for aluminum, E = 70.4 GPa and v = 0347, can be 



TABLE AL Peak loads and (aading/unfoadinA rates iked m the area 
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load range 
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High 
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High 
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30 


Low 


It 1-120 


at 


10 


Low 



1562 



J. Mater, Res.. Vol T, No. 8. Jun 1902 



cx? m ~i t TbcrcooTrTC ni 



M J kt.ut /*r* m Ktr\r> 



** I2'30bd TdJLOi ** 

JUN22 2005 14:06 FR CISTI ICIST 



TO 18027690501 



P. 22 



W.C Oliver and GuM. Wwr An improved technique for determining hardness ami 
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used to eliminate the unknown value of £ r in Eq. (18), 
This results in quicker convergence. The values of the 
elastic constants we use for the diamond indenter are 
Ei - 1141 GPa and *, = Ol07. The psoceduie also pro- 
duces A vs Ac data which ate used in the determination 
of the area function at large depths. 

The problem with using ahnaingm to extend die 
area function to small depths is that because of its tow 
hardness, small indentations in aluminum require very 
smaU loads, and a limit is set by die luce resolution of 
the indentation system. This problem can be avoided by 
making the small indentations in fused quartz, a much 
harder, isotropic material available in optically finished 
plate form from the General Electric Corporation- The 
standard procedure we use for determining the area 
function involves making a series of indentations In 
fused quanz using all the peak loads shown in Table AL 
Hie contact areas and contact depths arc then determined 



using Eqs. (IS) and (8) in conjunction with the reduced 
modulus computed from die elastic constants for fused 
quartz, £ = 72 GPa and i/ s 0.170. Because the ma- 
chine compliance is known from the aluminum analysis, 
there is no need for iteration. 

lb finally establish the area taction, the A vs ^ data 
obtained for both the al uminum and the fused quartz are 
& as a group to Eq. (19). The area function at large 
depths is then determined primarily by the aluminum 
data while that at small depths rs determined by the fused 
quartz data. A convenient fitting routine is that contained 
in the Kakidagraph™ software for Apple Macintosh™ 
co m p ut ers, A weighted procedure can be used to assure 
that data points with small and large magnitudes are of 
equal importance. For the loads outlined in Table AI, 
rite mtntmum contact depth is about 15 nro and the 
maximum about 4700 nm, so the area function is good 
in the range m between. 
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Abumcl— Mcchaaical Stress ill antcrconncctiaQ* tS A arobtcitt of 
growing Importance fal VLSI djwicaa. Open circuits from metal cracV- 
ina; and voiding awd ttart circuits from hillock* pne atrtts-celBtcd pax* 
nemeaa. The o*-*gt«« of thU stress an discussed Including intrinsic 
stresses from the sjotnests of the films sad thtrmalty Induced stresses. 
A measurement technique based on tfce determination at water cur- 
vature with a (a«er mauning device U utilised to directly OlCftnrt tbc 
film Areas in titu as a function of temperature during thermal cycling, 
Ttw chuoges in stress observed during thermal cycle* are interpreted 
quantitatively and foccfeamsat* that lead to plastic datamation' and 
tfaehr relationship to hillocks art discussed- l» tfc* stress vert us tenv 
pcrarore irajnsurements, several regions have bees identifitd uattudine 
elastic nod plastic behavior both uitaercamm*w1on and (castas, (be 
yield strength, recrystailUaUon, gain growth, hardema^, and solW- 
stste reaction*. The effects of deposition conditions an these regions 
are also esamlned. 



I. Introduction 

SEVERAL failure mechanisms associated with me- 
chanical stress are becoming increasingly important in 
thin films used for tnterconnecuona in integrated circuits. 
Open circuits from metal cracking and voiding, short cir- 
cuits due to hillocks, and both open and $hon failures due 
to electromigraiion are phenomena related to mechanical 
stress. Temperature is important because the conductors 
and dielectrics have relatively large thermal mismatches 
with silicon resulting in stress. Solid-state reactions and 
changes in structure normally involving volume changes 
may also occur causing variations in stress during pro- 
cessing [1], [2]. in addition, many deposition processes, 
particularly the low-tcmpcratuxe processes, result in films 
with high intrinsic stress (3] because of defects or struc- 
tural mismatch between the firm and substrate [4]. The 
origin of these stresses has been reviewed by Hoffman 
[SI Campbell [6], [7J. and Kinosita [8], 

The stresses introduced by these causes are often large 
and may exceed the strength of the film; the result can be 
cracking or, if an interface tails, dclaminauon. The rets* 
live softness of aluminum easily results in deformation. 
Aluminum interconnections are typically covered wttn a 
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dielectric layer deposited at elevated temperatures to pro- 
tect them from corrosion and to enable multiple levels of 
wiring. This dielectric expands and contracts at a different 
rate than the aluminum doting thermal cycles, which 
causes change in the stresses. In addition, the dielectric 
can have its own large stress that can affect the encapsu- 
lated aluminum [91. Thermal cycling also occurs during 
the standard anneal or alloy in H 2 + N 2 at 4S0*C toward 
the completion of microchip fabrication. The stresses 
generated can lead to the formation of protrusions or hil- 
locks that adversely affect the integrity of the dielectric 
layer. A review of the mechanical properties of thin films 
with emphasis on theory and strain relaxation was pub- 
lished, including experimental data from lead films to il- 
lustrate the theoretical results f 10], 

Studies of stress have helped to achieve a better under- 
standing of hillocks, and because hillock suppression is a 
primary advantage of layered films (1 1H13J, this knowl- 
edge of internal stress has assisted in explaining the im- 
provements [\\ .An additional benefit was the detection 
of change in thin-film stress when film composition and 
structure changed (1), f 14]. This has proven to be a very 
useful advantage because the precise temperature at which 
processes such as compound formation or rccry stalliza- 
tion occurs can be determined. 

II. Stress and Hillocks 

A generalized stress-strain curve is illustrated in Fig. 1 
tor comparison to actual measurements. The elastic re- 
gion is linear, reversible, and time independent. When 
stress becomes high enough, the curve is no longer linear 
because the deformations arc becoming permanent and, 
when force is removed, the resulting sample length is not 
the same. This can be seen in the shift alter unloading in 
Fig. 1(a). When the applied stress creates a permanent 
strain of 0.2 percent in soft materials such as aluminum, 
the yield strength has been reached; this strength or stress 
is not always the same in the measurements thai follow 
because it is dependent on deposition conditions and com- 
position. 

An understanding of the resole of applying force to thin 
films is not easily obtained, Dtffusional flow at the grain 
boundaries was thought to be the predominant mechanism 
for stress relaxation f 16} in evaporated films while only 
under tension, but similar behavior was not observed in 
sputtered films. Another Study pointed out that grain- 
boundary diffusion alone cannot totally relax the stress in 
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Fig. I. Centralized Btrcswitrein curves for w» Mvmxrairo «ih>y (15). (*) 
Two ttile regions, elastic and ptutic. (b) Yield «rcnjth defiaed as a 
0.2 peiceat deviation, 

a film attached to a substrate {17}. The dominant mecha- 
nism in thin films of aluminum Is probably thermally ac- 
tivated dislocation glide {19] even above naif the melting 
point- Diffusion*! now does play a role because disloca- 
tion motion will occur at lower stress when dislocations 
climb in addition to glide; such a climb occurs through 
the thermally activated diffusion of vacancies. 

Hillocks appearing in thin films arc a result of com- 
pressive plastic deformation. If a technique is used that 
reduces hillock density, it will not generally eliminate 
plastic deformation; chin films will always behave etasti- 
cally and plastically when force* ate applied. The mor- 
phology or structure of the film can be changed to 
strengthen it, thereby minimizing dislocation motion, 
which would reduce hillocks if dislocations play a role m 
their generation. A better approach utilized in mulrilay- 
ered films is to create a surface that would allow Uniform 
plastic deformation but block the development of hillocks 
at any localized point. The material roust be strong enough 
to preclude deformation bat does not need to be ductile 
because the strain rates ate low. Strength is needed to pre- 
vent material from breaking or .diffusing through the sur- 
face layer and to distribute the stresses evenly across the 
aluminum film* There will stifl be plastic deformation, 
but it would occur more uniformly. 

Several materials can be used to achieve a strong sur- 
face layer. One is the deposition of a high strength ma* 
. terial free of pinholes, such as titanium (111. A second 
could be a thick aluminum oxide grown on the film; alu- 
minum has the advantage of a relatively hard amorphous 
oxide as can be seen in nanoindentation measurements 
[19]. This approach has been investigated (20}-[22] and 
is being refined [23]. To eliminate hillocks, however, the 
native oxide must be 0.4 *m for a 0-7-nm initial alumi- 
num film thickness [23]. The third is an Ion implant into 
the aluminum surface to create a highly damaged layer 
that strengthens the film f24H2S] because of the many 
dislocations and also possibly because the ions them- 
selves produce point-detect Hardening. When dislocation 
density is high, the individual dislocations interfere with 
each other, which prevents plastic deformation and effec- 
tively strengthens the film. This is the same mechanism 
that strengthens a material when it is work hardened. The 
disadvantages of this technique are that the ion-implant 



damage can anneal out after a high-temperature cycle and 
the dectrotnigration lifetime at room temperature can be 
lower. 

til. Origins of Stress in Thin Films 

Mechanical stress in thin films can occur during the 
deposition process (mtrinsic) or is generated from changes 
in temperature. That from thermal cycling is either the 
result of differences in the thermal expansion of the ma- 
terials or stnicxuxal changes within the materials and can 
be quite large. Because the difference to the linear expan- 
sion coefficients is large for aluminum versus silicon and 
me strength of aluminum is relatively small, film defor- 
mation b difficult to avoid. The low degree of strength 
can actually be aft advantage because the stresses that lead 
to anamination can be prevented from developing be- 
cause plastic deformation will relieve stress. This is a 
common problem with thick films of refractory metals be- 
cause the severe intrinsic and thermal stresses result in 
strong forces at the interfaces leading to delamination. An 
additional advantage of aluminum is its ability to reduce 
S1O2 and form its own oxide, which further improves 
adhesion* thereby facilitating the use of thicker films. 

A. Intrinsic Stress 

Intrinsic Or initial stress from the synthesis of the film 
can be as nigh as several gigapascals in sputter-deposited 
films. The initial stress may be the result of differential 
thermal expan s io n caused by changing temperatures dur- 
ing deposition or may be a result of the growth of the film 
(26]. Without measurements of the substrate temperature 
as a function of time during deposition, it is difficult to 
account for rite unrntemrnnal variations in temperature; 
therefore, in the context of this study , intrinsic stress will 
be the initial stress if no substrate hearing is used. 

Stress that develops in films as they am formed can be 
caused by different mechanisms. Changes in volume dur- 
ing the formation of the film can lead to intrinsic stress;' 
for example, when SK> 2 is grown on the silicon surface, 
the change in volume causes large compressive stresses 
that are not easily relieved m low deposition temperatures 
{27]. Intrinsic stress can also develop from the coalescing 
of islands of material during the initial stages of film for- 
mation {26}; stress was found to exist in these islands be- 
cause of lattice mismatches and, when they grow to- 
gether, the stress changes. 

In sputtering, intrinsic stress is dependent on argon-gas 
pressure and was modeled by an atomic peening mecha- 
nism [28] based on the elastic rebound of neutralized 
sputtering tons (argon ions) from the target. These ener- 
getic gas atoms bombard the deposited film as it conden- 
ses on the substrate [29], 130] and these reflections are 
dependent on such parameters as die atomic-mass num- 
ber. The larger elements have a greater cross-sectional 
area that results in a higher number of sputtering atoms 
(such as argon) reflecting from die target and subse- 
quently striking the film being deposited. 
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a Thermally induced StrvsS 

When the temperature of a thin film on a substrate is 
changed, a strain is introduced because of the variation in 
thermal expansion 

«ntm = \ {acta - a^) dT (l) 

J To 

where a^ m and arc the linear-expansion coefficients 
of the film and substrate, respectively. It is possible to 
obtain a polynomial expression for a from the literature 
[31], [321 and then to integrate. Over limited temperature 
ranges and at high temperatures, mis expression can be 
approximated as 

eiita 88 (<*am ~ AT (2) 
which assumes that the linear-expansion coefficients are 
independent of temperature. 

Stress in the thin film can be determined directly from 
the experimental measurement by obtaining the substrate 
curvature. A calculation of the strain can be used to then 
predict the stress in thin nlmS- 

In mululayered films, plastic deformation is possible in 
one film while the others behave elasticity The elastic 
component can be predicted; however, plastic behavior 
wilt make it difficult to model the measured stress* Point* 
by-point subtractions of the curvature can be made, how- 
ever, if a separate measurement is made for each individ- 
ual layer. This has its limitations in that the individual 
films could oxidize or react with the substrate, therefore 
behaving differently than the multilayer film. For combi- 
nations of conductors and insulators such as aluminum 
and silicon dioxide, this technique has worked well * 

Changes in stress from thermal heating can also be 
caused by structural changes in the material. For exam- 
ple, different crystal structures can have different vol- 
umes, and tilms can soften from recrystallizaiion or be- 
come stronger from work hardening. If such a change 
occurs during thermal cycling, it will be reflected in the 
stress behavior as a function of temperature. 

IV. Measuring Film Stress 
The stress in thin (Urns can be measured by optical in- 
tcrfcromctry [33], [34], X-ray diffraction to measure the 
lattice spacing (35] or to measure the substrate curvature 
[36H39], canttlevered bending beams [40], and laser 
scanning [18)» [41 }. [42). The most accurate and simplest 
is to measure the substrate curvature using the laser- scan- 
ning method. This is a derivative measurement in that the 
measured change in displacement of the reflected laser 
beam is proportional to the change in the slope of the tan- 
gent line to the surface of the wafer. This change in dis- 
placement can also be thought as being proportional to the 
variation in the angle between the incident laser beam and 
wafer surface. Accuracy is approximately 5 percent be- 
cause calibration is achieved with standard minors. Pre- 
cision, however, is 2 x 10~* m -1 or an equivalent radius 
of curvature of 5 km, which translates to a stress of 1.5 



MPs In a l-*im film on a 500-jtm substrate. Because this 
ts approximately 1/100 the yield strength of aluminum, 
the measurement will be accurate enough to observe, for 
example, clastic/ plastic strain transitions. 

Stress was measured in all films after deposition and as 
a function of temperature by means of a laser to determine 
in situ during thermal cycling the change in the radius of 
curvature. The apparatus is detailed in the literature (1], 
[lft|, [43]. Stress in thin tilms is calculated by measuring 
the change in wafer curvature and using the following 
expression {44], (45]: 




where a«m is the stress in the thin film, £ wi> and are 
the modulus of elasticity and Poissoo's ratio for the sub- 
strate, r„ b and r« ra are the substrate and film thickness, 
# a is the initial radius of curvature before film deposition, 
and R is the measured radius of curvature. 

Note that this calculation requires a change in the radius 
of curvature, the elastic properties of the substrate, and 
the thickness of the substrate and film. It is not necessary 
to know the mechanical properties of the film, assuming 
it Is of uniform thickness, which indicates that, indepen- 
dent of composition and niicrostructure. stress can stilt be 
obtained. The resulting measured stress in the thin alms 
is then analyzed and modeled based on any known me- 
chanical properties. Stresses in the substrate and thin film 
are further Retailed in the literature [46], [47]. 

Aluminum-based metal films were deposited onto < 100) 
oriented 10-20 uSl • cm p-typc silicon wafers by mag- 
netron sputtering at room temperature. A thermal oxide 
100 run thick was grown prior to thin-film deposition. This 
SiOj layer was grown to prevent epitaxial regrcwth of the 
aluminum on silicon and to prevent dissolution Of the sil- 
icon substrate in the aluminum during thermal processing. 

Stress in films of Al-$i and Al-Cu was measured, and 
the results are plotted in Fig. 2. The temperature was 
raised to 450° C to simulate the effects of semiconductor 
processing. Measurements were obtained every 12.5 de- 
grees. This proved advantageous because several details 
previously not observed (43] became apparent. The ordcr- 
of-magnitudc improvement in accuracy and the control- 
lable heating rates also contributed to better measure- 
ments. The heating rate was 5*C/min and the ambient 
atmosphere was air. 

The initial as-deposited Stresses in the two films in Fig. 
2 differ partly because die degree of thermal conduction 
from the substrate to the holder resulted in variations in 
the deposition tempecatu res. These intrinsic stresses can 
also be affected by the temperature of the source* which 
can vary as a result of differing powers, presputtering 
times, and cooling water* If the intrinsic stress in the film 
docs not reach the yield strength during cooling down 
from the deposition temperature, the point at which stress 
becomes zero is dose to the actual deposition temperature 
assuming that very tittle stress is generated during film 
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formation. Such mechanisms as ioa bombardment Mid the 
latent heat of condensation increase die effective deposi- 
tion temperature; bombardment can also change the plas- 
tic properties of the film* For depositions above approxi- 
mately I50 Q C, some plastic flow in the aluminum films 
after cooling down from the deposition will keep the 
stress-free point near 150*C as has been observed by oth- 
ers 149]. 150]. 

The Sims in Fig. 2 were deposited at room temperature. 
The temperature at which stress in the aluminum films 
changed from compressive to tensile, however, is be- 
tween 80° and 100° C, which implies that the average 
substrate temperature daring deposition may have been 
90° C. This heating is primarily from the kinetic energy 
that the aluminum obtains when sputtered by the Ar + ions 
from the target. These ions have energies of several elec- 
tronvolts [51]. An earlier study reported that the temper- 
ature of deposition for spattered aluminum on unhealed 
substrates increased by approximately 90° C 152], which 
agrees well with the data. Stress in evaporated aluminum 
deposited at room temperature was also determined by 
measuring the curvature of die wafer and was found to be 
near zero. This correlates well with the fact that evapo- 
rated aluminum atoms have typical energies of 0.2 to 0.3 
cV [51], Evaporated films of copper, titanium, and chro- 
mium onto cold substrates were found to have an ratrmsic 
stress 1401, but this is probably because the temperature 
of the surface needs to be higher for the atoms to become 
mobile. 

To explain the many changes in slope in the following 
temperature versus stress curves, the plot Tor aluminum 
(with 1-percent silicon) is separated into regions and sim- 
plified. The bullc moduli of elasticity of aluminum, silicon 
{531, and titanium are listed in Table I as are their melting 
points and hardness for comparison, and their use for dun 
films is appropriate 145], The bard^s/strengm «f the ti- 
tanium aluminice and its relatively high melting point are 
indicative of the good success in strengthening the alu- 
minum imerconncctrons- It will be difficult tor hillocks to 
form in layered Al /Ti films as a result of the hardness of 
the atuminide because plastic deformation in the alumi- 



num will be forced to be uniform. Other alumimdes hav- 
ing a large number of aluminum atoms (such as WAl u ) 
are probably not as strong. 

V. Regions in the Tbwpbraturb Versus Stress 
Curves 

There are many stress variations in thin films of alu- 
minum as a function of temperature. Categorizing the dif- 
ferent sections of stress behavior led to the identification 
of at least eight regions. With each variation in slope, 
there is a change In the mtcrostxucturc or film behavior. 
For example i during plastic deformation, a deviation to- 
ward zero stress indicates that the film may have become 
softer. With a more rapid move away from zero stress, 
the film may have become stronger if the rate of change 
does not exceed the elastic Stress expected from the ther- 
mal mismatch between the film and substrate. Fig. 3 plots 
stress as a function of temperature and was chosen be- 
cause seven of the eight regions are apparent. This num* 
her of regions in one film is not common. The data points 
were removed for clarity, and the remaining line is a least 
squares fit of the data points. The numberfi denote the re- 
gions. Region 3 is not shown because it was not observed 
in pure films of aluminum but was seen in nearly every 
alloy studied and is due to solid-state reactions [1], (2], 
[»]. 

A. Elastic Region 

The first region is the basic elastic behavior observed 
in most materials. Stress was found to be reversible in 
that, by heating and cooling between 25° and 100°C, the 
path followed was linear and would retrace itself, which 
is typical of elastic behavior. Tins stress can be calculated 
via Hooke's lav* and the strain resulting from the differ- 
ence in the expansion coefficients. Combining Hooke's 
law with (2) results in 



(cw-ot^AT (4) 



where is film stress and E/O - *) is the biaxial 
modulus of elasticity in the thin film. This expression 
models stress resulting from changes in temperature. Be- 
cause there can be an initial intrinsic stress, however, it 
would not correctly determine the absolute stress, and a 
better expression would be the derivative of stress as a 
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fig. 3, Regions in chc 



function of temperature as 



Mm _ f ^ ] 

r ' L(i - •.)! 



(5) 



where Aa is the difference in the expansion coefficient 
between the thin film and the substrate. Based on the val- 
ues Of £/{ 1 — a) for (lti) aluminum and « tor (ill) 
aluminum and (100) silicon from Table 1, a theoretical 
slope of -2.27 MPa/*C correlates relatively well with 
the measured value. An additional advantage of this mea- 
surement is that the difference in the expansion coeffi- 
cients or the modulus of elasticity for thin films can be 
determined knowing the other 

B. Recrystallizatiun and Grain Growth 

The property changes that follow after thermally cy- 
cling materials that have been cold formed or deformed 
occur as three mutually distinguishable processes— recov- 
ery, recrystaUization, and grain growth (54]* Recovery 
results in some dislocation rearrangement and annihila- 
tion, but the original grain boundaries are basically re- 
tained. RecrystaUization begins with the formation of new 
unstrained dislocation-free grains that grow in the strained 
matrix or deformed material; its primary driving force is 
the stored energy from plastic deformation but depends 
on time, percetiuge deformation, temperature of defor- 
mation, purity, and grain site £55]. Grain growth occur- 
ring in completely recrystailized metals is a continued 
growth of recrystailized grains. These three processes re- 
quire time, but the time is exponentially dependent on 
temperature. When the temperature is sufficiently high, 
these processes appear to be instantaneous. 

Region 2 in Fig. 3 represents a large reduction in strain 
corresponding to a softening of the film. This softening is 
the result of recrystaUization and grain growth and has the 
effect of removing damage caused by the sputtering pro- 
cess. RecrystaUization eliminates dislocations from the 
crystallites of aluminum us can be observed in the samples 
in Fig. 4. Dislocation networks can be seen in the TEM 
before but not after thermal cycling to Similar be- 




Pig. 4. Rec*>«*Uiz*ion of oluimmun. (*) Aa deposited, (b) After hearing 
to200*C. 



havior has been reported but was attributed not to rccry s- 
talitzarioo but only to grain growth [56]; however, rc- 
crystaliizaiton in bulk aluminum induced by small 
amounts of plastic strain has been observed 157]. Both 
processes produce larger grains. Normally, recrystaUiza- 
tion is followed by grain growth which absorbs and drives 
the dislocations in die matrix to the grain boundaries. Be- 
cause dislocations strengthen a film, the newly formed 
grains have reduced strength and hardness. 

After the dislocation- tree crystallites form, grain growth 
will continue and further soften the films because there 
are fewer grain boundaries co stop dislocation movement 
within tiie crystallites. The difference between rccrystal- 
lizauon and grain growth is that recrystaUization occurs 
when the material is strained or deformed, whereas grain 
growth relies on grain-boundary energy. Growth can oc- 
cur if the grain boundaries do not have urtpurtties (such 
as oxides), which typically pin the boundaries after the 
initial grains form during deposition. 

Stress is frequently observed to decrease after some 
plastic deformation as in the Al-Cu curve in Fig. 2. Re- 



BEST AVAILABLE COPY 



on * j 



TThrrcooTrrc ni 



GAftDNEK AND PUHNi STaCSS AS A FUNCTION OF T£MP£H4TU*F W M RLMS 



crystallization caused by plastic deformation * assumed 
to occur in region 2 based on two additional observaiioas. 
Pint, strengthening Ac material with titanium foyers 
raised the temperature at which die stress began jo drop. 
Grain growth driven by grain-boundary energy should not 
be suppressed by titanium. Second, when die films were 
cooled ra liquid nitrogen < which has flte effect of placing 
the films under tension and reducing plastic flow) and then 
warmed hack io room temperature, the mm will be at the 
compressive yield strength at a lower temperature. As a 
result, plastic deformation will begin at a lower temper- 
ature. This procedure reduced the temperature at which 
«*ftcning of the film began by 10" to 15°C, which would 
not occur if there was not a dependency on plastic defer- 
mat job . 

Grain growth can be seen in the TEM in Fig. 4, The 
average grain size as deposited is approximately 0.2 pm, 
which increases to 1 pm after heating to 200*C When 
the size is small compared to film thickness, growth is 
driven by the reduction of gram-boundary area and en- 
ergy v which results in a columnar-grain structure wherein 
all grain boundaries extend from the top to the bottom of 
the film and all grain-houndary planes axe roughly per- 
pendicolar to the plane of the film. The resulting mean 
grain size tends to be comparable to film thickness. 

Abnormal or secondary grain growth is also possible 
and was observed in films of a sandwich-layered At-€b- 
Ai structure (5Sf ; the size of the resulting grains was often 
as large as 100 #mvlt was reported that, if precipitates of 
CuAla arc present, such grain growth Is not possible, and 
a graph of the temperature at which this growth occm* 
was observed *o follow the solid aolttbOity of copper in 
aluminum. Large grains were also found in films of Al- 
90-percent Cu and were correlated with stress {59}. The 
grains were not abnormally large when the films were de- 
posited on alumina diaphragms. A study of a sandwich- 
layered Al/Mg/Al reported abnormal grain growth (iSOJ 
and work with an Al/Cr/Al structure also reported 50 to 
100 jim grains (611; however, the effects of stress have 
not been studied and hillocks were a problem. Stress plays 
an important role in the generation of these grains (SSJ. 
The growth of large aluminum grains also occurred in 
films deposited at temperatures above 56Q*C £62}. 

C Elastic/ Plastic Tnittsiiioa 

Region 4 in Fig. 3 exhibits a gradual change from the 
slope reflecting elastic behavior to a slope of zero where 
stress does not increase although the rising temperature 
results in strata caused by thermal mismatch. This occurs 
because the film is soft and* in some aluminum films de- 
posited under ultrahigh vacuum conditions, mis transi- 
tion region was not observed because of the high purity 
of the films; instead, stress remained constant and low 
immediately after reerystallizatiou resulting from plastic 
deformation. Fig. 5 is a measurement of the stress in such 
a film mat was synthesized using a high deposition rate 
and a background pressure in file lower 10~* mbar range 
obtained by baking the samples and pumping the sputter- 




Kg. 3. temperature verag arcs* of a Upta Ai-Sl tampfe fabricated using 
* better vacuum *nd Mgb deposition rate*, Tbe background pressure dur- 
ing depntition wa« J * 10~* rater and tfca depoiittofl jwwor «ro »t the 
10 kW raaiiroura resetting » * spmxettag time of* ram. The totowh 
rmr of the tttteteuef «#t 12 rpm. 

tog system over a period of two days. Using a better vac- 
uum results in softer aluminum films because trace 
amounts of oxygen or water in the deposition system react 
with aluminum. This results in precipitates of A1 2 Q, or in 
small amounts of A1 2 0, at the grain boundaries. The lack 
of precipitates reduces die aluminum yield strength as can 
be seen in the first hearing in Fig, 5- The second heating 
has a transition region at 200 *C, but k Is small. The alu- 
minum is seen to have a higher yield strength, but this is 
because the film does not recrystallize during the second 
thermal cycle and, if the sample were to remain at a tern- 
perauuc above 200°C, the stress would relax. 

The yield strength or stress of a material is dependent 
on many factors including impurities and temperature and, 
as a result, it is difficult CO obtain reliable data for com- 
parison. Table JJ lists values for high- and low-purity alu- 
minum. Yield strength is defined as the stress required to 
produce a specified small permanent strain (normally 0.2 
percent) as illustrated in Fig. I, and such a strain can oc- 
cur with the appropriate change in temperature. Equation 
(2) can be used to calculate the change in temperature that 
corresponds to this amount of strain and, based on the 
di if e recces in the expansion coefficients for aluminum and 
silicon. tOG°C was obtained for a 0-2 -percent strain. To 
determine the point where the yield strength has been ob- 
tained, the temperature at which the stress would be zero 
as if no recrystallizatioa occurred is then required plus the 
change in temperature corresponding to the clastic strain. 
If rectystaliizatxon has occurred, zero stress would be at 
approximately 150*C This translates to u line similar to 
that in Fig. 1(b) representing elastic behavior drawn from 
25Q*C (because 100°C corresponds to 0.2 percent) with 
a Slope equivalent to the clastic region ( -2 XI MPa/ °C ) . 
The yield strength of the film being measured, therefore, 
is the intersection between the 0.2-pcrcent strain line and 
the curve and this occurs at 300°C; therefore, the fourth 
region appears between 200* and 300* C. 

In soft materials such as aluminum, the boundary be- 
tween the elastic and plastic regions is not well defined, 
and a slight permanent deformation may take place even 
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for very small stresses [15], This region, therefore, is re- 
ferred to as an elastic/plastic transition region, and its ap- 
pearance after recrystaUizauon depends on the strength 
and purity of the thin film. Such a region is always present 
during subsequent thermal cycles because of the film tran- 
sition* from eliuttic to plastic behavior as it changes from 
under tension (positive stress) to under compression (neg- 
ative stress). 

D. Compressive Yield Strength 

The measured stress in the film no longer increases with 
temperature in region 5 of Fig. 3, end the film is plasti- 
cally deforming. Elastic strain is constant and the yield 
strength of aluminum has been exceeded. This strain con- 
tains an elastic and a plastic component as 

*m« = t> + t* (6) 

When stress is lower than the yield stress, strain is purely 
elastic and, when h is higher, there win also be plastic 
strain that will hold the total film stress constant. 

Because the film is under compression, this is the pe- 
riod during which hillocks arc most likely to develop* es- 
pecially as a result of the high temperature , Hillocks were 
observed to slowly grow just above the recrystatltzauon 
temperature (150*C), which is closer to being catego- 
rized as region 4*, therefore, region 5 is not the only sec- 
tion where compressive stresses are high enough to form 
such protrusions. The plots in Fig. 6 are the reflected in- 
tensity of the laser beam as a function of temperature of 
the Al-Si sample used in Fig. 5. This intensity is seen to 
dmp- above 250*C partly because of the formation of hill- 
ocks on the aluminum surface. Because plastic deforma- 
tion will always occur when this stress level is obtained, 
hillocks can be eliminated most efficiently by controlling 
the deformation in such a manner that it becomes uni- 
form. This can be accomplished with layered films. 

Typically . the strength of materials begins to decrease 
at approximately half the melting point of the material. 
This is known as thermal softening and, in aluminum, the 
temperature around which it begins is 200°C. A 5 a result, 
yield strength should become lower with rising tempera- 
ture as illustrated in Fig. 3. The dashed line denotes the 
predicted yield strength as a function of temperature and, 
if the film had come under compression sooner, it would 
be expected to follow this curve. In very soft films and at 
higher temperatures, stress can drop very close to zero. 
In addition, if the temperature is held constant at 450°C, 
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stress can be seen to decrease exponentially over a period 
of hoots; its final value depends on film composition and 
structure. These phenomena are discussed further in our 
previous work [in. 

£ Elastic Behavior During Cooling 

When temperature is lowered, compressive stress is ini- 
tially reduced dasticaily (region 6) because there is an 
elastic and a plastic component of strain in the film as 
shown in (6). The film will eventually pass through a tran- 
sition to tensile stress and will behave elasttcaUy until 
yield strength is again reached which, at this elevated 
temperature, will be lower than when it was under 
compression. The elastic to plastic transition is more 
gradual because the strength of the aluminum wilt be in- 
creasing. The slope of the clastic region will also be 
slightly different because the elastic moduli will be smaller 
for silicon and especially for aluminum at these temper- 
atures. 

r. Tensile Yield Strength 

The yield strength of aluminum under tension is reached 
In region 7. Its absolute value is higher than the strength 
under compression partly because hillocks that have 
formed do not easily draw tack into the film. Plastic de- 
formation lakes place as the temperature is further re- 
duced. This can also occur immediately following depo- 
sition if the deposition temperature is high enough for ten- 
sile stress to exceed yield strength after the samples are 
cooled. 

This region of stress was not observed in earlier mea- 
suemeats 156] and can be used to determine the yield 
strength at room temperature (19]. Plastic behavior still 
occurred in the earlier measurements because the final 
stress obtained after cooling to room temperature was 
lower than would be expected if behavior was purely clas- 
tic; however, based on the higher strength, the films may 
have been hardened. One possible mechanism that can 
strengthen a film is precipitation hardening, and such pre- 
cipitates can be formed by trace amounts of oxygen or 
water. The background pressure in the earlier investiga- 
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tion [56] was 1 x 1<T* torr wuh a deposition raw of 360 
nm/min. This base pressure results in approximately one 
monolayer per second of gas impurities deposited on the 
samples or 3 percent of die film, which is high enough for 
AlaOa precipitates to form and affect the film. 

Another effect noted in the earlier work £561 and in this 
investigation [1], [19] is a higher yield strength with re- 
duced film thickness partly as a result of hardening that 
occurs with smaller grain sizes. The film thickness is typ- 
ically related to the grain size and smaller grains wfll re- 
sult in a stronger 61m because dislocations do not propa- 
gate as easily through the matrix, In addition, mere can 
be a true thickness effect independent of grain size that 
could be studied by etching back a thick film after grain 
growth. The stress measurement of two spntter-deposited 
thicknesses of aluminum clearly reveal a difference in Fig. 
7. Stress in thin aims of lead were also found to have a 
very strong film-thickness dependence (351. If the final 
value in the stress versus temperature measurement, which 
typically represents yield strength, is plotted as a function 
of thickness, thinner films of aluminum are seen to be 
stronger, in addition, nanoindcutafcon of these films dis- 
plays a similar trend in strength versus thickness as plot- 
ted in Fig, S. 

It has been proposed that stress relief normally occurs 
through dislocation slip where the limiting step is the gen- 
eration of dislocation loops [56], and the stress required 
to generate these dislocations is assumed to be an inverse 
function of grain size. In our recent work on me stress 
that develops in thin films of tungsten £19], however, the 
strength of the film was observed to dimmish with de- 
creasing film thickness, which contradicts the earlier pro- 
posal- Another alternative [63] is that the relative elastic 
modulus of the substrate'as compared to the thin film plays 
a role in driving the dislocations way or toward rhe film* 
substrate interface that can be modeled by the imaging of 
dislocations; however, a good quantitative model is not 
available to predict the effect on strength. Models in which 
die dislocations are pinned at the surface have been de- 
veloped to help explain the dependency of stress on film 
thickness [45]. 

G. film Hardening 

For region 8, several explanations have been proposed 
concerning the cause of the film hardening, but it is un- 
clear which is the most valid. One is that precipitation 
hardening can strengthen the film as illustrated in die tem- 
perature versus stress measurements m Fig. 2 for Al-Cu. 
As the film cools, the solid solubility of copper in alu- 
minum decreases (0 T 045 atomic percent at 225~C) and, 
at a certain level, excess copper will form the theta-phase 
CuAt 3 . This phase will result in precipitates that effec- 
tively strengthen the film. The temperature at which this 
becomes apparent is 17S°C, which is lower than typical 
temperatures for accelerated clcctromigration tests. This 
may indicate that die accelerated tests are being con- 
ducted on metallization when die copper is in solid solu- 
tion, whereas the operating conditions are at t e mpe r at ures 
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where precipitates have formed. Stress measurements of 
Al-Cu have also been recently measured by Townsend, 
and similar hardening is apparent in some of hie results 
[501, !<>*J- Unfortunately, no measurements of aluminum 
without copper were obtained for comparison. 

In the At-Si films, silicon will precipitate out as the 
tcrmjeratnre drops, and mis can result in limited harden- 
ing- The difference is that CuAl 2 tends to form numerous 
small precipitates while silicon can produce nodules as 
large as 1 jtm is size partially because of the very fast 
diffusion rate of silicon in aluminum. 

Because hardening is also observed, although to a. lesser 
extent, in films that do not have copper or silicon (Fig. 
7), precipitation hardening cannot be the only reason. 
Dislocations interacting with one another may have a more 
profound effect at lower temperatures. Other effects due 
to dislocations may also contribute [1). In addition, when 
copper or silicon in solid solution forms CuAl 2 or silicon 
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precipitates, a change in volume can increase tensile stress 
if there is a reduction in volume. 

VI. Conclusion 
When mechanical stress is applied to materials, elastic 
and plastic deformation occurs. Elastic behavior can be 
modeled quantitatively whereas there are various mecha- 
nisms thai can lead to plastic deformation and are less 
easily modeled. In aluminum aims, yield strength is eas- 
ily reached because of the large difference m the thermal 
expansion coefficient between it and silicon and because 
of its relative softness and low melting point. As a result, 
plastic deformation involving dislocation motion occurs. 
The manufacturmj! process itself can affect the behavior 
of the films to terms of deposition temperature* number 
of passes in from of each target {13], and the quality of 
the film. Smooth surfaces free of hillocks can be devel- 
oped using additional layers of materials such as titanium 
and should be taken into consideration when designing a 
film. 

The stress in dun aluminum films and its alloys was 
measured in situ during thermal cycling la a laser-based 
stress apparatus as a function of - temperature and then 
modeled. Both intrinsic and thermally induced stresses 
were observed and examined. These stresses are respon- 
sible for the formation of hillocks and for such problems 
as thin-film delaminatlon, cracking, and voiding. During 
thermal cycling, at least eight regions of stress were iden- 
tified, including elastic behavior, recrystallization, plastic 
behavior, yield strength, film hardening from precipi- 
tates, grain growth, and compound formation from solid- 
state reactions. A more complete and comprehensive study 
of several other aluminum alloys are presented in our other 
studies II]. J2J. [141. 
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